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By balancing both accuracy and efficiency, density func- 
tional theory (DFT) (IJ 13 and its time-dependent exten- 
sion (TDDFT) f3^, ^41 have become arguably the most em- 
ployed methods in computational physics and chemistry. 
DFT and TDDFT are based on rigorous theorems, which 
reformulate many-electron quantum mechanics using the 
simple one-electron density as the basic variable of inter- 
est rather than the complicated many-electron wavefunc- 
tion. In this letter we prove that the theorems of TDDFT 
can be applied to a class of qubit Hamiltonians that are 
universal for quantum computation. In a similar spirit to 
DFT and TDDFT for electronic Hamiltonians, the theo- 
rems of TDDFT applied to universal Hamiltonians allow 
us to think of single-qubit expectation values as the basic 
variables in quantum computation and information the- 
ory, rather than the wavefunction. From a practical stand- 
point this also opens the possibility of approximating ob- 
servables of interest in quantum computations directly in 
terms of single-qubit quantities (i.e. as density function- 
als). Additionally, we demonstrate that TDDFT provides 
an exact prescription for simulating universal Hamiltoni- 
ans with other universal Hamiltonians that have different, 
and possibly easier-to-realize two-qubit interactions. 

We begin by briefly reviewing TDDFT for a system of N- 
electrons described by the Hamiltonian 

H{t) = l^f + l^w{\ri-fj\)+fvir,t)n{r)d\, (1) 

.^^ zm .^j J 

where p^ and are respectively the position and momen- 
tum operators of the ith electron, w(|f / — fj|) is the electron- 
electron repulsion and v(r,r) is a time-dependent one-body 
scalar potential which includes the potential due to nuclear 
charges as well as any external fields. n{r) = Y4 5{r — Vf) is 
the electron density operator, whose expectation value yields 
the one-electron probability density. The first basic theorem of 
TDDFT, known as the "Runge-Gross (RG) theorem" |4|, es- 
tablishes a one-to-one mapping between the expectation value 
of n{v) and the scalar potential v(r,^) and therefore through 
the time-dependent Schrodinger equation, a one-to-one map- 
ping between the density and the wavefunction. The RG the- 
orem implies the remarkable fact that in principle, the one- 
electron density contains the same information as the many- 
electron wavefunction. The second basic TDDFT theorem 
known as the "van Leeuwen (VL) theorem" |5 1 gives a pre- 
scription for constructing an auxiliary system with a different 
and possibly simpler electron-electron repulsion w'(|f; — 
which simulates the density evolution of the original Hamil- 
tonian in Eq.[T] When w'(|f / — rj|) =0, this auxiliary system 
is referred to as the "Kohn-Sham system" |2| and due to it's 



simplicity and accuracy, is in practice used in most DFT and 
TDDFT calculations . 

It is not obvious that the RG and VL theorems extend to 
qubits, which are distinguishable spin 1/2 particles. We now 
prove analogous RG and VL theorems for a system of N 
qubits described by the universal 2-local Hamiltonian (61 [T], 

i= 1 i= 1 i= 1 

(2) 

Here, of, of, of are Pauli operators for the ith qubit, hi{t) 
are local applied fields arbitrarily chosen along the z-axis and 

and J^ij^i are two-qubit interaction terms respectively 
parallel and perpendicular to the direction of the fields. The 
above Hamiltonian describes an open chain of N qubits ar- 
ranged in a one-dimensional array, with each qubit interact- 
ing with its nearest neighbors. In Refs. OITl, it was shown 
that by appropriately tuning the local fields in Eq. [2j one can 
realize any two qubit gate, which in turn can be employed 
to perform universal quantum computation when combined 
with singe-qubit rotations. Although Eq. [2] describes quan- 
tum computing with a fixed two-qubit interaction, by applying 
appropriate local fields we can generate effective Hamiltoni- 
ans describing systems with tunable two-qubit interactions as 
well. This can be viewed as a time-dependent version of the 
widely employed method of Gadgets {VT\. In Eq.|2] the case 
where /^^^ = 7^"-^^ yields the Heisenberg Hamiltonian which 
describes exchange coupled spins in solid state arrays or quan- 
tum dots in heterostructures |8|. The situation /^^^ 7^ ^l+i 
yields the XXZ Hamiltonian, used to model electronic quoits 
on liquid Helium |9| or solid-state systems with anisotropy 

due to spin-orbit coupling llQl . while the limit 7^"-^^ = 
yields the XY model describing superconducting Josephson 
junction qubits ifTB . 

We now state the equivalent RG theorem for quantum com- 
putation with the Hamiltonian in Eq. [2j 

Theorem - For a given initial state |v^(0)) evolving to \ w{l)) 
under the Hamiltonian in Eq. and with and Jj^i^i 

fixed, there exists a one-to-one mapping between the set of 
expectation values {cr^, crl, ...cr^} and the set of local fields 
{/zi,/z2, .../^iv} over a given interval [0 ,t]. 

Here, we have defined af = ( V^(0 1 of \ \\f{t)) as the expectation 
value of the component of the ith qubit along the field direc- 
tion (z-axis). A detailed proof together with a more rigorous 
discussion of the conditions on the theorem are provided in 
the supplementary material. The RG theorem implies that the 
set of local fields can be written as unique functionals of the 
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Figure 1 : Runge-Gross theorem for a 3 qubit example - The set of 

expectation values {dp cjI, ...cr^}, defined by the the Bloch vector 
components of each qubit along the z-axis in (a), is uniquely mapped 
onto the set of local fields {/zi,/z2, ...%} in (b) through the RG the- 
orem. Then, through the Schrodinger equation, the set of fields is 
uniquely mapped onto the wavefunction. These two mappings to- 
gether imply that the N-qubit wavefunction in (c) is in fact a unique 
functional of the set of expectation values {Cp cri, ...cr^}. 



set of expectation values {(7^ O"!, ...cr^}, as illustrated in the 
first part of Figure [T] Since the solution to the time-dependent 
Schrodinger equation is unique and /^^^ and jf^^i are fixed, 
the wavefunction is a unique functional of the local fields, i.e. 
\\l/{t)) = |i/a[/zi,/z2, .../?iv](0)' where the square brackets de- 
note that i/A is a functional of the set {/zi,/z2, .../?iv} over the 
interval [0,t]. This fact, combined with the RG theorem al- 
lows us to state a corollary, which is the first central result of 
this letter: 

Corollary - There exists a one-to-one mapping between the set 
of expectation values {(7^ O"!, ...cr^} and the N-qubit wave- 
function \^f(t)) on the interval [0 ,t]. 

The above corollary implies the counterintuitive fact that the 
full N-qubit wavefunction, which lives in a 2^ dimensional 
Hilbert space, is a unique functional of only the N compo- 
nents of each qubit along the z-axis. i.e. 



(3) 



Although the RG theorem does not tell us an explicit func- 
tional form for i/A, it has profound conceptual implications 
from a quantum information perspective. At first glance, it 
might appear that the set {(7^ O"!, •••O'^} contains much less 
information than the full wavefunction, since projective mea- 
surements needed to obtain {cTp al, ...CT^} would typically 
imply that information about non-commuting obesrvables is 
lost. However, since the wavefunction completely specifies 
all properties of the system, Eq.|3]implies that even properties 
depending on non-commuting observables such as entangle- 



ment and phase information are in fact unique functional of 
the set of expectation values {cr^ O"!, ...cr^}. 

From a practical standpoint, the RG theorem implies that 
all observables can directly be approximated as function- 
als of single-qubit expectation values, without regard for 
the wavefunction. Although the set of expectation values 
{o"p cri, ...O"^} in principle contains all of the quantum infor- 
mation in i/A, extracting this information directly is not always 
straightforward, although in some cases it is. As an example, 
consider a computation involving only one flipped qubit rela- 
tive to the other N —\ qubits having an opposite orientation. 
An explicit entanglement functional (as measured by concur- 
rence 1 13 1) between any two qubits labeled k and / can be 
written very simply as 



Eki[olGl,...Glj]{t)- 



1 



n 

m=kJ 



N 



(4) 

Interestingly, this particular entanglement functional is time- 
local, since it depends only on the set {cr^ a"!, ...cr^} at a given 
instant in time and so we may write Eki[G^,G2^ ...G^]{t) = 
Eki[Gl{t)^G2{t)j...G^{t)]. In the more general case, ob- 
servables may be non-local in time and depend on the set 
{g^,G2^ ...G^} over an entire interval [0,t]. Although the func- 
tional in Eq.|5]is time-local, it is spatially non-local, since the 
entanglement between qubits k and 1 depends on the compo- 
nents of all of the other N — 2 qubits. If one considers two 
flipped qubits instead of one, the entanglement functional be- 
comes complicated and non-local in both space and time due 
to dependence on phases in the wavefunction (see supplemen- 
tal material). Understanding the spatial and temporal non- 
locality of density functionals in electronic structure theory 
is a very active research topic (Ml [151, and naturally arises 
here in TDDFT for quantum computation as well. 

We now turn to the second fundamental theorem of TDDFT 
for universal computation, a VL-like theorem for qubits: 

Theorem - Consider a given set of spin components 
{(7^, (jI, ...O"^} obtained from the wavefunction \ y/{t)) evolved 
under the Hamiltonian in Eq. |2] There exists ( see supplemen- 
tary material for certain conditions) a Hamiltonian with dif- 
ferent two-qubit interactions denoted J^f^i and /-'i^i and dif- 
ferent local fields {h[ , /z2, -"h^i^}, which evolves a possibly dif- 
ferent initial state \ V^'(O)) to a different final state | V^'(O) 
that the condition {g[^,G2,...g^^} = {a^, a"!, ...cr^} is satis- 
fied on the interval [0,t]. 

Here, we have defined g'^ = {\ir' (t)\Gf\\if' (t)) . The VL theo- 
rem allows us to obtain the set {cr^, O"!, ...cr^} by simulating 
the evolution with an auxiliary Hamiltonian having different 
two-qubit interactions and hence a different (and possibly sim- 
pler) wavefunction evolution as illustrated in Figure [2] This 
opens the possibility of simplifying computations by con- 
structing simple approximations to the auxiliary fields as func- 
tionals of single-qubit expectation values, in the same sense 
that the exchange-correlation potential of electronic DFT and 
TDDFT is approximated as a functional of the one-body den- 
sity in the Kohn-Sham scheme. 
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Figure 2: Van Leeuwen theorem for a 3 qubit example - The set 
{(7p (jI, ...CJ^} (a) obtained from evolution under Eq.[2] is uniquely 
mapped to a new set of fields {h[,h2, .-.h'j^} (b) for a Hamilto- 
nian with different two-qubit interactions. Evolution under this new 
Hamiltonian returns the same expectation values {dp (72, ...CJ^}, al- 
though the wavefunction is different and hence projections of the 
Bloch vectors along other axes are in general different (c). 



The proof of the VL theorem also gives a mathemati- 
cal procedure (see supplementary material) for engineering 
the exact auxiliary fields {h\^h'2^...h'j^} which reproduce a 
given set {cr^, O"!, ...cr^} under a different two-qubit interac- 
tion. As a simple demonstration, we use this procedure to 
numerically simulate a 3 -qubit Heisenberg Hamiltonian us- 
ing an XY Hamiltonian as the auxiliary system (Figure [3]). 
For the simulation, the system is prepared in the initial state 
|i^(0)) = -ij(|011) + |101) + |110)), where |1) and |0) are 

eigenstates of with eigenvalues -1 and 1 respectively. In 
the Heisenberg Hamiltonian, /^^^ = J^"-^^ = Ji^i+i and we 
choose J\2 = J23 = 0.5. We apply a pulse of the form 
h{t) = 0.61^^1 (-l)"+isin [(2n - l)t] to the first qubit and 
hslt) = 0.61:^^1 (-l)^"sin[2nt] to the third qubit. The time- 
dependent Schrodinger equation is solved numerically and the 
set {(7^,(72,(73} is read out during the evolution. Details of 
the simulation are provided in the supplementary material. 

and we choose 



For the auxiliary XY Hamiltonian, /■ 
different and non-uniform couplings in which =1.2 and 
= — 1. Using the VL theorem, we engineer the auxil- 
iary local fields {h[^h2jh^2} which using a non-uniform XY 
interaction, reproduce the set {(7^, (7^, O"!} obtained from the 
original evolution under the uniform Hesienberg Hamiltonian. 
As seen in Figure [3] the auxiliary local fields are quite differ- 
ent from the original local fields applied to the Heisenberg 
model, but simulate the set of components {(7^, al, O"!} cor- 
rectly, i.e. {cr(^, (72^, cr3^} = {(7^, al, O"!}. In the language of 
electronic DFT, the XY model in our simulation is analogous 
to the "Kohn-Sham system" and the set {h\^h'2-,h'^} play the 
role of the exact Kohn-Sham potential as a density functional. 



The VL theorem could also be useful from an experimental 
perspective, as it allows one to engineer different pulses which 
perform the same computations, but using different two-qubit 
couplings. 

The RG and VL theorems place TDDFT for universal quan- 
tum computation on a firm theoretical footing and open sev- 
eral exciting research avenues. The development of approxi- 
mate density functionals has been essential for the success of 
electronic DFT and TDDFT and will be in quantum compu- 
tation and information theory as well. Functionals that per- 
form actual computational tasks or subroutines such as the 
quantum Fourier transform are immediate research directions. 
The collapse of the computational complexity class hierarchy 
is of course not expected, and therefore finding functionals 
that carry out complex quantum computational tasks is ex- 
tremely unlikely. Nevertheless, understanding how TDDFT 
functionals can approximate quantum algorithms is an open 
direction. Density functionals for strongly correlated lattice 
and spin systems have been recently proposed 1 16 - 19] and 
could be applied to several problems of relevance in quantum 
computing. In Refs. |[T6lfT9l local density (LDA) and gener- 
alized gradient approximations (GGA) for one dimensional 
Hubbard chains and spin chains were derived from exact 
Bethe ansatz solutions and could readily be applied to solid- 
state quantum computing or perfect state transfer protocols 
in spin networks |20|. Functionals can also be parametrized 
from numerical simulations of one-dimensional qubit systems 
using time-dependent density matrix renormalization group 
methods (TDMRG) [21], in an analogous fashion as quan- 
tum Monte Carlo simulations of the uniform electron gas have 
proven invaluable in electronic DFT 1221 . 

Another important research direction will be the general- 
ization of DFT and TDDFT to other universal Hamiltonians 
and models of quantum computation. For instance, Ref. i23l 
discussed the use of TDDFT for obtaining gaps in adiabatic 
quantum computation. In Ref. 1291 , groundstate DFT was 
used to study relationships between entanglement and quan- 
tum phase transitions, while Ref. |31| explored DFT from a 
complexity theory perspective. We are also exploring an ex- 
tension of the TDDFT theorems to models with directly tun- 
able two-qubit interactions, as opposed to the fixed interac- 
tions discussed in this letter. 

Useful discussions with S. Mostame, J. D. Whitfield, S. 
Boxio, M. H. Yung and J. Parkhill are greatfully acknowl- 
edged. We thank NSF award PHY-0835713 for financial sup- 
port. 



Supplemental Material 

A. Proof of the Runge- Gross theorem 

In this section we will first consider a proof of the RG the- 
orem for Hamiltonians of the form, 

H{t) = + of of) + I4af + thim^, (5) 

i<j i<j i=\ 
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Figure 3: Simulating the Heisenberg Hamiltonain with the XY Hamiltonian - Pulses of the form hi{t) = 0.6£^^^(— 1)"+ sin [(2n— l)t] 
and h^it) = 0.6^,^=1 (— l)^"sin [2nt] are respectively applied to the first and third qubits of a uniform Heisenberg Hamiltonian (a). The time- 
dependent Schrodinger equation is then solved exactly numerically and the evolution of the set {dp (72, (73} is read out in (b). The VL theorem 
gives us a prescription for constructing different auxiliary fields (c), which simulate the evolution of the set {0^,02,0^} correctly as seen in 
(d), but using a non-uniform XY interaction instead. (Time is measured in units of ^) 



5 



which reduces to Eq.[2]of the main text in the Hmit of a one- 
dimensional array with nearest-neighbor coupHngs and open 
boundary conditions. We will see that it is possible to formu- 
late the RG theorem of time-dependent current density func- 
tional theory (TDCDFT) for the more general class of Hamil- 
tonians in Eq. [5] but for TDDFT one must stay with the more 
restricted form in Eq.|2] 

The proof begins with the equation of motion for the expec- 
tation value of the ith qubit along the field direction (z-axis), 



(6) 



where (O) = {\\f{t)\d\\\f{t)) denotes the expectation value of 
an arbitrary operator O at time t and | \ir{t) ) is the wavefunction 
evolved on an interval [0,t] from a given initial state |v^(0)), 
under the Hamiltonian in Eq. [5] Development of the commu- 
tator in Eq.[6]yields, 



dt 



^ Gf=2tjtim^j)-{^i^f)) 



(7) 



where we have introduced the Pauli raising and lowering op- 
erators = 6"^ ± in the second equality. Defining 



Jki 



1 



-2JtM^J-alaf) = -JtAKK-^k^n (8) 



as the operator that generates the "current" of flowing from 
the ith qubit to the kth qubit, Eq. [7] takes the form of a local 
conservation law, 



dt 



cyf = -l,{jki)^ 



(9) 



k^i 



This arises from the fact that the Hamiltonian in Eq. [5] con- 
serves the total component of all N-qubits along the field di- 
rection, i.e. it is readily verified that. 



(10) 



This is analogous to the situation in electronic structure the- 
ory, where the local continuity equation 



dt 



n(r,f) = -V-j(r,f) 



implies a global conservation of particle number 



d_ 



J n{r,t)d' 



r=^N = 0, 



(11) 



(12) 



where N is the number of electrons in the system. 
We now consider a "primed" Hamiltonian 

i<j i<j i=\ 

(13) 



which has the same two-qubit interaction terms as the Hamil- 
tonian in Eq.ls] but a different set of local fields {h\^h'2^ ...h'j^}. 
Let Iv^'(O) denote the wavefunction evolved from the same 
initial state | V^(0)), but under the Hamiltonian in Eq. 



13 



The 

equation of motion for the expectation value of the ith qubit 
along the z-axis under this primed Hamiltonian is 



dt 



N 

-Y^^jkiYj 

k^i 



(14) 



where we define {6)' = {Y {t)^6\\if' (t)) as the expectation 
value of an arbitrary operator O with respect to the primed 
wavefunction. 

In what follows, we assume that the local fields in both 
the primed and unprimed systems are equal to their Tay- 
lor expansions within a finite radius of convergence around 



0. 



i.e. hi(t) = l-Aj,_^hi(t) 



t' and simi- 



r=o 



larly h[{t) = LJ^q 



t=0 



(the assumption of Tay- 



lor expandability is discussed below). We will now pro- 
ceed to show that if the set of fields {h\^h2-, ...hjq} dif- 
fer from the set of fields {h'^^h'2^...h'j^} by more than 
a global field which is the same for all N qubits, 
the set of currents {{jn) , ...{ji?) , ...{jn-\,n)} and 
{{jnY, {juY, ...(As)', ••'{jn-i.n)'} win necessarily be differ- 
ent. The condition that the two sets of fields differ by more 
than a global field, is equivalent to the statement that there 
exists a smallest integer m ^ such that the set 

)^—(h,{t)-h\m,=o,^{h2{t)-h'2mt=o, 

(%(0-/!^(0)|f=o}7^{C}, (15) 



■■■ dt 



where here {C} is a constant set of N elements that are all 
the same. i.e. the Taylor coefficients of the local fields in the 
primed and unprimed systems will differ at some order. 

Next, we write down the equation of motion for the dif- 
ference of the currents between the ith and kth qubits in the 
primed and unprimed systems: 



dt 



aJki) - ikY) = i{[H{t)Jki])-i{[H{t),U)' . (16) 



Since both systems evolve from a common initial state | Va(0) ) , 
we have at ^ = 0, 



dt 



i{iki)-{ikiy)\t=o=i{vm [{m-H\o)),j\i] \xifm 



= 4(vA(0)|4-|vA(0))(A/i,-(0)-A/z^(0)). 



(17) 



Here, we have defined Ahi{t) = hi{t) — h[{t) as the difference 
between the unprimed and primed fields acting on the ith qubit 
and similarly, Ahk{t) = hk{t) — h'j^{t) . 

% ^ 4{dldf + aldj) = ^f{d+dr + 6^ at) (18) 
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is similar to a local kinetic energy operator, describing the 
total transfer of cr^ between the ith and kth qubits. From 
see that if the condition in Eq. JTS] is satis- 



Eq. 17 



we 



fiedform = 0, the sets {(712), (jb), ...(723), ...(7>-i,iv)} and 
{(/12)', 0"i3)', •••{h?>)'i '•'{Jn-i.n)'} will become different in- 
stantaneously later than ^ = (with a restriction on the vanish- 
ing of (i/A(0)|f^; |i//(0)) discussed below). If the condition in 
Eq. [15] instead holds for some m > 0, we differentiate Eq.[T6 
m times to obtain, 



{{jki) - {jki)')\t=0 



= 4(vA(0)|4-|vA(0))^(A/zKO-A/i;^(0)|.=o. (19) 

From here we see that if the set of local fields eventually differ 
at any order, the set of currents must as well. This establishes 
the RG theorem of TDCDFT: For a fixed initial state |i/a(0)), 
there is a one to one mapping between the set of local fields 
and the set of currents, up to a globally constant field. 
We now discuss the three main conditions of the theorem: 

1) The expectation values (i/A(0)|f^; |i/A(0)) must be non-zero 
for at least one pair of qubits k and i, whose local field dif- 
ferences Ahi{t) and Ahk{t) are different for at least one instant 
on the interval [0,t] . This is a fairly mild restriction on the set 
of admissible initial states, |v^(0)). For instance, consider a 
worst case scenario, in which all the fields {/ii,/z2, .../^iv} and 
{/Zp/z2, ..Mn} differ by a constant field, except for hi and h[ 
which differ by a different amount from the others at only one 
instant in time on the interval [0,t]. In this worst case, the re- 
striction means that (va(0)|7i/| ^/(O)) must be non-zero for at 
least one value of i, where / = 1,2, ...,A/^. In the more general 
case, where the sets {h\^h2-, ...hj^} and {h[^h2^ ...h^} differ 
for several qubits or on finite time intervals, this restriction is 
even less severe. 

2) The elements of the sets {/zi,/z2, .../?iv} and {h[,h2^ ...h^} 
must be analytic functions of time. i.e. equal to their Taylor 
expansions within a finite radius of convergence. In quantum 
computing, this is not a very severe restriction, as one typi- 
cally constructs pulses which are well behaved functions. This 
restriction does not even exclude sudden switching, which is 
the case when applying idealized pulses to perform single- 
qubit rotations. 

3) The theorem establishes a one to one mapping between the 
set of currents and the set of local fields up to a globally con- 
stant field, we will denote C{t), which is the same for all N 
qubits. If one applies periodic boundary conditions, the extra 
symmetry fixes the value of C(^) and the mapping is one to 
one between the fields and currents with no ambiguity. For 
open boundary conditions, C{t) remains arbitrary, which cor- 
responds to an arbitrary term C{t)Y4 af = C{t)Gf^^^i in the 
Hamiltonian. If one begins in an initial state which is an 
eigenstate of O"^^^^^/, this term is simply a c-number and adds a 
trivial global phase to the wavef unction. This is typically the 
case when one begins in a computational basis state. How- 
ever, if one starts in a superposition of states with different 



values of of^^^^, the term C(t) of^^^^ yields a nontrivial coher- 
ence between these states. Such coherences would be measur- 
able for observables with non-zero matrix elements between 
states of different o^^^^i- These observables would therefore 
not be uniquely determined by the current when considering 
open boundary conditions. 

We now turn to the RG theorem of TDDFT, 
which is discussed in the main text. From Eq. |9] 
we see that it is possible for two different sets 
of^ currents ^{{jn),{jn),---{h?,),-{jN-\,N)} and 
{{jnY Ah^y i'-{h^)' i-'{Jn-\,n)'} to correspond to the 
same set of spin components {cr^, O"!, ...cr^}, if there exists a 
set of current differences 

{5jl2, 5ji3, ...5j23, ---^jN-lM 

^ {{{jn)-{jn)')Mh)-{h)'). 

... {{h)-{h)').A{h-l,N)-{jN-l^N)')}. (20) 



such that. 



k^i 



(21) 



for some i. For a one-dimensional chain with open boundary 
conditions and nearest-neighbor couplings as in Eq. |2] such a 
set never exists, as illustrated in Figure]?] Thus, for this case, 
no two sets of currents can yield the same set {cr^, a"!, ...cr^}, 
and through the RG theorem of TDCDFT, no two sets of fields 
{/zi , /z2 , . . ./?iv} and {h\^h'2-, •••h'^} differing by more than a con- 
stant can yield the same set {cr^, cri, ...O"^}. This establishes 
the RG theorem of TDDFT for the Hamiltonian in Eq. |2] of 
the text. For more general geometries, such as in Figure |5] 
it is possible to find two different sets of currents such that 
Ylk^i ^ jik = 0. For these geometries, the RG theorem of TD- 
CDFT holds, however that of TDDFT does not. From the con- 
tinuity equation of electronic structure (Eq. [TT]), we see that 
we can add an arbitrary transverse vector field 5 j (such that 
V • 57 = 0) to the electronic current, without altering the value 
of §jn{r,t). A set {5712, Sjb, ...5j23, ...5jV-i,iv} satisfying 
the condition in Eq.|2T]is analogous to a purely transverse cur- 
rent. 



B. Proof of the van Leeuwen theorem 

The proof of the VL theorem begins with the equation of 
motion for the current under the evolution of the Hamiltonian 
in Eq.|5] 



j^{jki)=i{[H{t)j\i]). 

Developing the commutator in Eq.[22]yields, 
d 



(22) 



dt 



{hi) = i^ki) + W + mdShikit), (23) 
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^2 



dt 



at ^3 



Jl2 



J23 



dt ^2 



^^3 



ii2 



^23 



at ^4 



J34 



at ^4 



^34 



Hamiltonian, 

i<j i<j i=\ 

(26) 

having different two-qubit interactions f-j- and Z^' and a dif- 
ferent set of local fields {h[^h2^ ...h^j^}. We allow the initial 
state I V^'(O)) to be different from the initial state | V^(0)), with 
the only constraint that the initial currents and the initial ex- 
pectation values of must be the same in the primed and 
unprimed systems, i.e. we only require that, 



{v\0)\jl\xif\0)) = {xifmj\MO)) 



(27) 



and 



Figure 4: Open Chain of 4 qubits with nearest-neighbor cou- 
plings - If we consider the set of currents shown in a), it will never be 
possible to find a new set of currents that will yield the same deriva- 
tives of (7^. For instance, as shown in b), we can find a new set of 
currents that will yield the same derivatives of on the middle two 
sites, but the derivatives at the ends of the chain will necessarily be 
different. 



where 5hik{t) = hi{t) — hk{t) is the difference between the lo- 
cal fields applied to the ith and kth qubits and we have defined 
the operators and as: 



af- 



mk^m 



(7/ 





y J^ a^ 















(24) 



and 







y /" a' 










y /" 











y /" 



m^k 



(25) 



^ki arises from the commutator of the current operator with 
the kinetic energy operator Y4<j Tij and is similar to the "stress 
tensor" operator of electronic TDDFT | 5 1. arises from the 
commutator of the current operator with the term jfjOfdj 
in the Hamiltonian and is analogous to the "internal force den- 
sity" operator due to the electron-electron repulsion in elec- 

repre- 



23 



tronic TDDFT. Both the terms (^,) and {^ki) in Eq.[: 
sent "internal forces" due to the two-qubit terms in the Hamil- 
tonian, while the term A{fj^i)5hik{t) represents an "external" 
driving force due to the applied local fields which couple to 
the one-qubit operators of. 

We now consider an auxiliary "primed" system, with the 



(V/'(0)|df|v/'(0)) = (v/(0)|df|v/(0)) 



(28) 



for all k and i. The equation of motion for the current under 
evolution of this Hamiltonian is, 



dt 



-A{%)'5h[,. (29) 



Here, the operators Z^-, f^'-, c^'- and ^ / • ar e defined exactly as 
for the unprimed system in Eq' s [8) [24| and [25]respectively, 



but with /-y and f-j replaced by J[j- and /, ' 



are 



We now assume that all quantities in Eq.'s |23] and [29 
equal to their Taylor series expansions with a finite radius of 
convergence around ^ = 0. Denoting the Taylor expansion of 
an arbitrary function f{t) as f{t) = Y4=of^^^ ^^d Taylor ex- 
panding both sides of Eq. [23j we find after equating coeffi- 
cients of equal powers of t. 



{m+ 1 = i^kir + {Air +41^ {nr-'Sh^. 

s=0 



(30) 



Similarly, for Eq. 29 we have 



\'m+l 



We now subtract Eq. |3T] from Eq. [30] 
the set of currents {(712), (/is), ...(723 



(31) 

and demand that 
S...(/iv-i,iv)} and 



{ (/12 ) (/i 3 (^23 (/a^- 1 ,iv )'} be the same in the primed 
and unprimed systems, i.e. we demand that the Taylor coeffi- 
cients (jki)^ and (j^.) ^ be the same for all m and for all qubit 
pairs k and i. This yields 



m— 1 



s=0 



s=0 



+ {^kir-{^i:r^{Air 



(32) 



for all k and i. We see that the left hand side of Eq. 32 
contains Taylor coefficients of 8h\j^{t) of order m, while the 
right hand side has only Taylor coefficients of 5h\j^{t) of or- 
der less than m and known quantities. Thus, when supple- 
mented with the condition in Eq. [27] Eq. [32] is a unique re- 



cursion relation for the Taylor coefficients of the local field 
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Figure 5: Open chain of 4 qubits with second nearest-neighbor couplings - If we add second nearest neighbor couplings as well (shown 
in a)), it will be possible to find two sets of currents that yield the same derivatives of o^. For instance, if one considers a new set of currents 
around the closed loop shown in b), such that {5 ju^S 713, 5 723} = constant, the derivatives of will remain unchanged. This is equivalent to 
adding a purely transverse term to the current. 
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differences 8h\^{t), which reproduce the given set of cur- 
rents { (712) , (/i3 ) , . • . (723 ) , . . . {jN-]ji) } using different two- 



qubit interactions f-j- and J-j . Eq. 



32 



can be used to construct 
the set local fields, {h[^h2j ...h^p^} up to an arbitrary constant 
field. For periodic boundary conditions, the arbitrary constant 
is again fixed by the extra symmetry and the fields themselves 
are uniquely determined. 

So far we have established a VL theorem for TDCDFT. 
In order to establish the VL theorem of TDDFT discussed 



in the text, we simply need to add the condition in Eq. [28 
With this additional constraint, it is clear from Eq's. |9] 



and 
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that if the constructed fields {h[^h2^ ...h^^} force 
the s^{(;;2)',(/;3)'--(/23)'--(/^v-i,iv)'} to be the same 
as {(/i2),(/i3),.-.(/23),...(/iv-i,iv)}, the sets {a[\a2\ ...o'^} 
and {(7^, a"!, •••cr^} must be the same as well. 

We now discuss two main conditions of the theorem: 

1) All of the quantities appearing in Eq.'s|23]and|29]as well as 
the sets {cr^, O"!, ...cr^} and {cr(^, '"(^n} must be equal to 
their Taylor expansions within a finite radius of convergence 
for the theorem to hold. This is a much more restrictive con- 
dition than for the RG theorem, which only requires that the 
sets of fields {/zi,/z2,.../?iv} and {/Zp/z2, .../^J^} be equal to their 
Taylor series expansions. This restriction arises in the VL the- 
orem of electronic TDDFT as well and approaches to circum- 
vent this condition have begun to be researched 1241 . 

2) For the entire set of fields {h[^h2^ ...h'^} to exist, (t^i)'^ 
must be non-vanishing for all pairs of qubits k and 1. This 
too is a more severe restriction on the class of admissible ini- 
tial states than in the R.G. theorem, which only required that 
(tki)^ be non- vanishing for certain values of k and 1. However, 
since we are free to choose |v^'(0)) so long as it satisfies the 
conditions in Eq.'s 27 and [28] we will often be able to choose 
an initial state such that {t^-) ^ 7^ 0. From a practical stand- 
point, we have also found in our numerical simulations that 
for vanishing (t^-) ^, we can add a small convergence factor 
to make the fields well behaved at the initial time with little er- 
ror in the overall propagation. This situation does not arise in 
electronic TDDFT for continuous systems, but similar prob- 
lems have been noticed when one defines electronic TDDFT 
for lattice systems ||25H28]| . Since the qubit Hamiltonians we 
consider in this letter are also discrete, it is not surprising that 
a similar situation arises. 



C. Entanglement as a functional of the set { dp , . . . cr^ } . 

In this section, we will discuss the construction of two-qubit 
entanglement as a functional of the single-qubit expectation 
values, {cr^, (7|, ...cr^}. We use the concurrence as a measure 
of the entanglement between any two qubits in an N-qubit sys- 
tem |[T3l . Since the concurrence depends on non-commuting 
two-qubit observables, we expect that it is in general very hard 
to construct exactly as a functional of the set {cr^, O"!, ...cr^}, 
which are simple expectation values of commuting single- 
qubit observables. We will see that this is indeed the case. 



but in the spirit of electronic TDDFT, one can hope to develop 
simple approximations. 

The two-qubit reduced density matrix (2RDM) for the kth 
and 1th qubits is obtained by tracing the full N-qubit density 
matrix over all other N — 2 qubits in the system. In this letter 
we consider only pure states, so the 2RDM is simply given by 



P)^/ =Tri,...k_i,k+i,...i_i,i+i,...N[|v^(t))(VA(t)|], 



(33) 



where Tr denotes a partial trace. Defining the "time-reversed" 
2RDM as 



Pki 



(34) 



the concurrence E^i is defined in terms of the eigenvalues A/ 
of the matrix pkipkl as. 



(35) 



Eki max(0, V A 



In Eq. 35 the eigenvalues Xi are arranged in decreasing order, 
i.e. Ai > A2 > A3 > A4. 

We first investigate the concurrence for a system which is 
restricted to the o^^^^i = ±(A/^ — 2) subspace. There is only 
one flipped qubit relative to the other N —I qubits. This is 
also known as the single-excitation manifold. We denote | /) = 
|00...010...00) as the computational basis state with the ith 
qubit in the state |1) and all other qubits in the state |0) (the 
O's and I's can be interchanged without changing any results). 
The N-qubit density matrix can be expanded in terms of the N 
computational basis functions as. 



\w{t)){w{t)\ = l,ant)aj{t)\j){i[ 



(36) 



From the above expression, we find the 2RDM for the kth and 
1th qubits to be 



Pkl 



( Zi^kAam^ 0\ 

ai{t)al 

al(t)ak \ak(t)\^ 

\ 0/ 



(37) 



In Eq. [37] p^i is expressed in the 2-qubit computational basis 
states, {|00),|01),|10),|11)}. From Eq.'s [34] and |35l we find 
the concurrence to be. 



Eki = 2\ai{t)\\ak{t)\. 



(38) 



In order to re-express in terms of {cr^, O"!, ...cr^}, we need 
to invert the matrix equation 



(39) 



where and a are column vectors formed from the sets 



{(7^, O"!, ...cr^} and {\a\ p , 1^2 1 , } respectively and M 

is a square matrix with the diagonal elements equal to -1 and 
all other elements equal to 1 . Carrying out the inversion and 
substituting the result into Eq. 38 yields Eq.|4]of the main text. 
Thus, we see that in the case of a single flipped qubit, it is very 
simple to obtain an exact entanglement functional. 
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As a more complicated example, we consider entanglement 
in the of^^^i = ±(A/^ — 4) subspace, which contains states with 
two flipped qubits. We denote = |00...010...010...00) as 
the computational basis state with the ith and jth qubits in the 
state |1) and all other qubits in the state |0). Here, we can 

expand the N-qubit density matrix in terms of these 
computational basis functions as, 

W{t)){xif{t)\= £ a*j{t)aki{t)\kl){ij\- (40) 

i< j,k<l 

Obtaining the 2RDM as before, we find the eigenvalues of 
pkipki to be 

Ai = { /( L k7(0P)( i: MOP) + l L ^..-(04(01}', 



(41) 
(42) 



and 



^4 = 1 /(I k7(0P)(I MOP)-l 

(43) 

The terms depending on sums over the coefficients' moduli 
squared, \aij{t)\'^, can be obtained fairly easily in terms of 
{(7^, (jI, ...ct^} using the expression 



cyf = l-l^\au{t)\'. 



(44) 



However, we see that Ai and A4 also contain the term 
I Y.i^k,l^ki{t)al-{t)\, which depends explicitly on phases in the 
wavef unction. The imaginary parts of the coherences can be 
obtained from the currents, which in turn can be obtained from 
time derivatives of {cr^, (7|, ...cr^}. However, the real parts of 
the coherences depend on expectation values of the kinetic 
energy operators, f^/, and in general will need to be approx- 
imated. Since the number of computational basis states in- 
creases with the number of flipped qubits, we expect the exact 
entanglement functional to become progressively more com- 
plicated as more qubits are flipped. This highlights the need 
for constructing simple approximate entanglement functional 
and will be explored in future work. 



D. Numerical propagation of the VL construction. 

In this section we discuss how the proof of the VL theorem 
can be used to numerically construct a set of auxiliary fields 
{h\^h'2^...h'^}, which reproduce a given set { CT^ , O"! , . . . O"^ } us- 
ing a different two-qubit interaction. This procedure was used 
to simulate a Heisenberg model using an XY model in the 
main text, and was demonstrated in Figure [3] and Figures |6] 
and |7] (below). 

In principle, Eq. [32] can be used as a recursion relation to 
construct the Taylor coefficients of {h\^h'2^ ...h'^} io arbitrary 



order, but in practice this proves to be numerically cumber- 
some. Instead, we use a formulation of the VL construction 
based on a non-linear Schrodinger equation. A similar con- 
struction was presented in ref. 1 24 1 for electronic TDDFT. 

We begin by numerically solving the time-dependent 
Schrodinger equation in the "unprimed system". 



dt 



\xif{t))=H{t)\xif{t)), 



(45) 



for a given initial state |v^(0)) and a given Hamiltonian H{t) 
which we wish to simulate (see below for simulation de- 
tails). From \yf{t)), we can calculate all relevant observ- 
ables, and in particular, we can calculate the set of currents 
{ ( / 12 ) , ( /i3 ) , . . . ( 723 ) , . . . {}n- 1 ,iv) } at each time-step. 

We then construct the set of fields {h[^h2^ ...h^j^} which re- 
produce this set of currents, but using a Hamiltonian with dif- 
ferent two-qubit interactions Jlj- and T^j' , 

i<j i<j i=l 

(46) 



This is done by numerically solving Eq. 29 for the 
set {h\^h'2-,---h'^} at each time-step, with the re- 
quirement that {(/i2),(/i3),...(/23),...(/iv-i,iv)} = 

{(/;2)'.(/;3)'.-(/23)'.-(/^-i,iv)'} for an k and 1. i.e. 
we solve 



dt 



(47) 



for an 5h\j^{t 
from the solution to Eq. 



where on the left hand side {jki) are known 
However, since k^ki)' 



45 



(f^'.)' depend on the auxiliary wavef unction \\\f^(t)), to solve 
Eq. |47]we must simultaneously solve the auxiliary system's 
time-dependent Schrodinger equation. 



d_ 
dt 



xif\t))=H\t)\xi/{t)), 



(48) 



where H\t) in turn depends on 5h'-j^{t). Thus, Eq. 48 is a 
non-linear Schrodinger equation and Eq.'s |47] and |48|repre- 
sent a set of coupled non-linear ordinary differential equations 
for the auxiliary wavefunction {y/^t)) and the set of fields 
{h\^h'2^...h'j^}. To solve this system of equations, we begin 
with an initial state |v^'(0)) which satisfies the conditions in 
Eq.'s[27]and[28l By enforcing Eq.|28] we ensure that our solu- 



tions to Eq.'s 47 and [48] will reproduce the set {cr^, O"!, ...cr^j 



in addition to tEe currents. From | V^'(O)), we can solve Eq. 47 
at ^ = 0, 



dt 



{jki)\t=,={^in=^^{^L^^^^^^ 



(49) 

for field differences 5/z.^(0). After making a choice for 
the arbitrary global field, we can construct the set of fields 
{/z;(0),/z^(0),.../z;^(0)} at r = and the Hamiltonian ^'(0). 
We then solve 



dt 



(50) 
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to obtain \ Y{dt)) at the next time- step. From \ Y{dt)), we 
obtain {h\{dt)^h'2{dt)^...h'j^{dt)} by solving 



dt 



{jki)\t=cit = {^^ty\t=cit + (^;,)i.=.. +4(f/,)'i,=,,54(Jo. 

(51) 

This procedure is continued at each time- step until we obtain 
{h[,h2, •••h'^} and \ ^f'(t)) on the entire interval [0,t]. 

For the simulation presented in the main text, H(t) is the 
Heisenberg Hamiltonian and H'(t) is an XY Hamiltonian 
with the chosen parameters. With the initial state |v^(0)) = 
^(|011) + |101) + |110)), we solve Eq. 45 using the fourth- 
order Runge-Kutta method. | ^f{t)) is propagated on a uniform 
grid with 10^ time-steps, each of duration dt = 1.5 x 10~^^. 
With \^f{t)), we can calculate derivatives of the currents be- 
tween all 3 qubits to be used in Eq. |47] We use the procedure 
outlined in Eq.'s 47- 51 to obtain \^'{t)) and {h\.,h'2^h'^} of 
the auxiliary XY Hamiltonian. For the auxiliary system's ini- 
tial state, we chose |i/a'(0)) = |i//(0)) = ^TjdOll) + |101) + 
1 110)), which satisfies the conditions in Eq.'s 27 and 28 since 



the initial currents vanish in both the primed and unprimed 
systems. We also fix the arbitrary global field by choosing 
h^it) = h2{t) for all t. Since | V^'(O)) is an eigenstate of 0^^^^^, 
this choice corresponds to trivially fixing the global phase 
of the auxiliary system's wavef unction and any other choice 
would yield identical expectation values of observables. 

As shown in Figure Is] the set {cr^, cr|, O"!} is faithfully re- 
produced hy \\if'(t)). In Figure |6] we show expectation val- 
ues of several other observables calculated with \ ^f{t)) in the 
left column and \ in the right column. Naturally, ob- 

servables that depend explicitly on the set {cr^, O"!, crU are 
the same in both cases, while those that do not will be dif- 
ferent. In particular, we see that the entanglement is the 
same in both cases, since both |v^(0) and |v^'(0) remain 
a superposition of states with one flipped qubit during the 
evolution. This means that the explicit entanglement func- 
tional in Eq.|4] holds, and since both wavefunctions produce 
the same set {cr^, cri, O"!}, they necessarily produce the same 
entanglement. As expected, the currents {(/12), (723)} and 
{(/12)'' (723)'} same for both wavefunctions, while the 

kinetic terms {(112), (^23)} and {(7^2)'' (^23)'} different. 



The same situation arises in electronic DFT, where the Kohn- 
Sham wavefunction reproduces the correct density and cur- 
rent, but the kinetic energy is in general different from that of 
the true correlated wavefunction. In Figure [7] we plot the ex- 
pansion coefficients of | \\f{t)) and | Va'(^)) in the computational 
basis {|011),|101),|110)}, which as expected are rather dif- 
ferent. It is also interesting to note that in our formalism, the 
operators for the current and kinetic energy are different in the 
original and auxiliary systems, since we let all two-qubit pa- 
rameters in the Hamiltonian differ. This situation is different 
than in electronic TDDFT, where the kinetic energy and cur- 
rent operators themselves are the same, although expectation 
values may be different in the case of the kinetic energy. This 
is important in our formalism, especially with regard to the 
current, since although {(/12), (723)} = {(712/' (/23)'}' 
finds that in general {(712), (723)} 7^ {(712)', (723)'}. 

E. Analogies between electronic TDDFT and TDDFT for 
quantum computation. 



Throughout the manuscript we have tried to stress the 
analogies between electronic TDDFT and TDDFT for sys- 
tems of qubits. The main relevant quantities in electronic 
TDDFT and the analogous quantities in TDDFT for quantum 
computation are summarized in Figure [8] The current and 
local kinetic energy for electronic TDDFT are written for a 
single electron, but the extension to N electrons is straightfor- 
ward, by integrating over N —\ coordinates. 

Despite the clear similarities, there are also important dif- 
ferences between qubit and electronic systems. Qubits are dis- 
tinguishable quantum particles and the wavefunction does not 
need to obey any particular permutational symmetry. In con- 
trast, electrons are indistinguishable fermions with a fully an- 
tisymmetric wavefunction. Through the Jordan- Wigner trans- 
formation |30|, a system of qubits can in fact be mapped into 
a system of spinless fermions, but we do not pursue this in the 
present work. Also, the pauli sigma operators obey a different 
commutator algebra than the electronic position and momen- 
tum operators. The resulting kinetic energy operator is a two- 
qubit quantity in TDDFT for quantum computation, while for 
electronic TDDFT it is a one-electron operator. 
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Figure 6: Observables in the Heisenberg Hamiltonian verses the XY Hamiltonian - a) The entanglement calculated directly with the 
wavefunctions \\j/{t)) (left column) and \ w\t)) (right column) is seen to be the same, as expected since the propagation is restricted to the 
subspace with one flipped qubit and Eq.Hholds. b) The currents {(712)5 (723)} (left column) and {(/l2)^ (723)'} (right column) are also the 
same, as a consequence of the VL construction, c) The kinetic energies {(712)5 (^23)} (left column) and {(7'/2)^ (^23)'} (right column) are 
very different, as is often the case in electronic TDDFT as well. 
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Figure 7: Wavefunctions in the Heisenberg Hamiltonian verses the XY Hamiltonian - The real parts (a) and the imaginary parts (b) of 
the expansion coefficients of = Y4<:j(^ij(t)\ij) (l^ft column) and = L^<J^^J•(0IO) (right column) in the computational basis. 

Although the wavefunctions are clearly different, they both reproduce the same set {dp crl, (73} throughout the evolution. 
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Figure 8: Analogies between electronic TDDFT and TDDFT for quantum computation - Relevant quantities in electronic TDDFT (left 
table) and the corresponding quantities in TDDFT for quantum computation (right table). 



